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Abstract 

An algebraic structure related to discrete zero curvature equations is established. 
It is used to give an approach for generating master symmetries of first degree for 
systems of discrete evolution equations and an answer to why there exist such master 
. symmetries. The key of the theory is to generate nonisospectral flows (A^ = A', I > 

I 0) from the discrete spectral problem associated with a given system of discrete 

evolution equations. Three examples are given. 

X) ■ 

I Introduction 

I ■ 

^ I The theory of integrable systems has various aspects, although the term "integrable" is 

0^ ' 
O , 

O , are one of those important aspects and have a deep mathematical and physical back- 

^ ! ground. When any special character, for example the Lax pair, hasn't been found for 



00 



a given system of continuous or discrete equations, among the most efficient ways is to 
consider its symmetries in order to obtain exact solutions. It is through symmetries that 
Russian scientists et al. developed some theories for testing the integrabiUty of systems of 
evolution equations, both continuous and discrete, and classified many types of systems 
O ■ of nonlinear equations which possess higher differential or differential- difference degree 
symmetries (for example, see [@). Usually an integrable system of equations is referred 
^ I as to a system possessing infinitely many symmetries P] @]. Moreover these symmetries 
form nice and interesting algebraic structures [Q. 

For a given system of evolution equations Uf = K{u), both continuous and discrete, a 
vector field (t{u) is called its symmetry if (j{u) satisfies its linearized system 

^ = KH i.e., ^ = [K, a] := K'[a] - a'[K], (1) 

where the prime means the Gateaux derivative. Starting from a Lie-point symmetry, we 
can often construct the corresponding explicit group-invariant solutions. A symmetry a 
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may, of course, depend explicitly on the evolution variable t. If a symmetry a of the system 
Ut = K{u) not depending explicitly on t is a polynomial in t, i.e.. 



(^{t,u) =^-pi{u), n>l, (2) 

then we have 

pi = [K,pi_i], l<i<n, (3) 

and 

(adi^)"+Vo = 0, where (ad;^)po = [K,po\- (4) 

Therefore the symmetry (0) is totally determined by a vector field po satisfying (^). This 
kind of vector field po has been discussed in considerable detail and is called a master 
symmetry of degree n of = K{u) by one of the authors (BF) in Ref. 

The appearance of first degree master symmetries gives a common character for inte- 
grable systems of continuous evolution equations both in 1 + 1 dimensions and in 2 + 1 
dimensions, for example, the KdV equation and the KP equation. The resulting symme- 
tries are sometimes called r-symmetries (for more information, see ^ for example) and 
usually constitute centerless Virasoro algebras together with time-independent symme- 
tries [0 PI 0- Moreover these r-symmetries may be generated by use of zero curvature 
equations or Lax equations ||10[ and the corresponding master symmetry flows may also 



be solved by the inverse scattering method |11| |12|. In the case of systems of discrete 
evolution equations, there exist some similar results. For example, many systems of dis- 
crete evolution equations have r-symmetries and centerless Virasoro symmetry algebras 



p!3| [jT4| [T^, and the inverse scattering method may still be applied in solving themselves 



and their master symmetry flows |T8[ |T^. So far, however, to the best of our 

knowledge, there hasn't been a systematic mathematical theory to explain why there exist 
r-symmetries for systems of discrete evolution equations and how we can construct those 
r-symmetries when they exist, from the point of discrete zero curvature equations. 

Throughout this paper, "master symmetries" is used to express the first degree master 
symmetries which generate r-symmetries. Our purpose is to give an algebraic explanation 
of the first question above and to provide a procedure to generate those master symme- 
tries for a given lattice hierarchy. The discrete zero curvature equation is our basic tool 
to give rise to our answer and procedure. The Volterra lattice hierarchy, the Toda lattice 
hierarchy and a sub-KP lattice hierarchy are chosen and analyzed as some illustrative ex- 
amples, which have one dependent variable, two dependent variables and three dependent 
variables, respectively. 

Let us now describe our notation. Assume that u = (ui, ■ ■ ■ , Uq)'^ where Ui = Ui(t, n), 1 < 
i < q, are real functions defined over R x Z (in the case of the complex function, the 
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discussion is similar), and let B denote all real functions P[u\ — P{t, n, u) which are C°°- 
differentiable with respect to t and n, and C°°-Gateaux differentiable with respect to u. 
We always write £■ as a shift operator and 

{E"^x){n) = x^"^\n) = x{m + n), where X : Z — >R, m, neZ. (5) 

Note that x^"^^ here doesn't mean the m-th derivative. Set B'^ = {(Pi, • • • , P^)^! -Pi G 
B, 1 < i < r}, and denote by all matrix operators $ = ($ij)rxr where the entries 
$y = ^ij{t,n,u) G B, and by V, all matrix operators depending on a parameter A: 
U = {Uij)rxr, where the entries Uij = Uij{t,n,u, X) e B for all A, being C°°-differentiable 
with respect to A. 

We will need a multiplication operator 

[n\:B^ B, P[u\ ^ [n]P[u], {[n]P[u]){m) = m{P[u]){m), (6) 

which is often involved in the construction of master symmetries. This avoids an unclear 
expression nP[u], which may also mean {nP[u]){'m) = n{P[v]){m). For example, it is 
absolutely clear that ([n]P[-u])(m) = mu{m — 1) + mu{m), when P[u] = E~^u + u. We 
also need a difference operator A = E — E~^, whose inverse operator may be defined by 

1—1 oo 

{A-\){n) = {{E - E-^)-\){n) ■= ^{ u{n + 1 + 2k) - J2u{n - 1 + 2k)), (7) 

k=—oo k=l 

where u is required to be rapidly vanishing at the infinity. Moreover we define 

(A-^a) = (l/2)a[n], i.e., (A"^a)(n) = (l/2)an, a = const. (8) 
Obviously we can find that 

{E - 1)-' = A-^(l + E-'), (1 - E-')-' = A-\E + 1), (9) 

and thus 

{E - l)~^a = a[n], (1 - E~^)~^a = a[n\, a = const., (10) 

which may also be viewed as a definition of two inverse operators {E—l)~^ and (l—E^^)^^. 
Note that here we have used the operator [n] so that two functions {E — l)^^a and (1 — 
E~^)~^a have the other clear expressions. The operators A~^, {E — 1)~^ and (1 — E^^)^^ 
often appears in the expressions of master symmetries and thus master symmetries are 
usually nonlocal vector fields belonging to B'^. 

In order to carefully analyze algebraic structures related to symmetries, we specify 
the definition of the Gateaux derivative X'[S] of any vector valued function X & B'^ at a 
direction S & B'^ as follows 



X{u + eS), (11) 

£=0 
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which imphes that X' is an operator from to S*^, and need the following two product 
operations 

[K, S] = K'[S] - S'[K], K,S e B\ (12) 
[f,g\{\) = f{X)9{X) - f{X)9'{X), f,9 G C°°(R), (13) 

where C°°(R) denotes the space of smooth functions defined over R. It is known that 
(09, [-,■]) and (C°°(R), [■,■]) are all Lie algebras. 

We now assume that t/ G V'' and the Gateaux derivative operator U' is injective 
throughout the paper. Let us consider the discrete spectral problem 

E(j) = U(j) = U(n,u,XU, 

(14) 

(f)t = V(f) = V{n, u, A)0, 

where V G V^'. Its adjoint system reads as 

E-^tjj = UiIj = U{n,u,X)ilj, 

^ljt = {EV)^lj = {EV{n,u,X))'4^- 
Their integrability conditions are given by the following discrete zero curvature equation 

Ut = {EV)U - UV. (15) 

If the operator equation ([TB|) is equivalent to a system of discrete evolution equations 
Ut = K{n,u), K G B'^, then it is called a discrete zero curvature representation of Ut = 
K{n,u). Evidently 

Ut = U'[ut] + f{X)U^, if Ai = /(A), 

where U\ = Therefore a system of discrete evolution equations Ut = K{n, u), Kg B'^, 
is the integrability condition of ([I^ with the evolution law Xt = /(A) if and only if 

U'[K] + fUx = {EV)U - UV. (16) 

Note that the injective property of U' is indispensable in deriving zero curvature rep- 
resentations of systems of evolution equations. The equation (|16D exposes an essential 
relation between a system of discrete evolution equations and its discrete zero curvature 
representation. It will play an important role in the context of our construction of master 
symmetries. 

The paper is divided into five sections. The next section will be devoted to a general 
algebraic structure related to discrete zero curvature equations. Then the third section 
will establish an approach for constructing master symmetries by the use of discrete zero 
curvature representations, along with an explanation of why there exist master symmetries 
for systems of discrete evolution equations. In the fourth section, we will go on to illustrate 
our approach by three concrete examples of lattice hierarchies. Finally, the fifth section 
provides a conclusion and some remarks. 
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II Basic algebraic structure 



We aim to discuss Lie algebraic structures of symmetries including master symmetries 
by using zero curvature equations. It is natural to ask what algebraic structure exists, 
related to zero curvature equations. To answer this question, we first plan to expose a Lie 
algebraic structure for the space V", C°°(R)). 

Let {K, V, /), {S, W, g) G {B^, V^ C~(R)), in other words, K, S are vector fields, V, W 
are r x r matrix operators and /, g are smooth functions. We introduce their product 

UK, V, /), {S, W, g)j = {[K, S], IV, Wj, [/, gj), (17) 

where [K, S], lf,g} are defined by (p^ , (|13D, respectively, and [V, Wj is defined by 

[V, Wj = V'[S] - W'[K] + [V, W] + gVx - fWx, (18) 

where [V, W] = VW — WV. The same product as (|TB|) has been introduced for the 
continuous case in EU . 



Theorem 1 (Lie algebra) The space {{B'^, V, C°°(R)), [■, ■]) is a Lie algebra, the product 
|-,-] being defined by jjiTlp, i.e., 

{{K, V, /), {S, W, g)j = {[K, S], IV, Wj, [/, gj) 

where 

[K, S] = K'[S]-S'[K], 

IV, Wj = V'[S] - W'[K] + [V, W] + gVx - fWx, 
lf,gjiX) = nX)g{X)-f{X)g'{X). 

The proof of the theorem will be given in Appendix A. Upon looking at the product 
a little bit more carefully, we can find that the Lie algebra ((S*^, V*", C°°(R)), |-,-]) has 
a Lie sub-algebra {{B'^,V^ ,0), [■,■]), for which everything corresponds to the isospectral 
case. Moreover the center of an element of this Lie sub-algebra is often Abelian. 

The above theorem exposes that a Lie algebraic structure hidden in the back of vector 
fields. Lax operators and spectral evolution laws. Usually we just touch Lie algebraic 
structures of vector fields while discussing symmetries. If we analyze symmetries from 
the point of zero curvature equations, it is natural that we need to find and handle Lie 
algebraic structure for all triples {K, V, f) G {B'^, V"", C°°(R)) where K, V and / are related 
to each other by zero curvature equations. In other words, we need to observe how two 
triples {K, V, /), {S, W,g) appeared in zero curvature equations connect with each other. 
The following theorem tells us that such a kind of connection can be reflected by the Lie 
algebraic operation of {B'^, V, C°°(R)) in Theorem |l|. Its proof can be found in Appendix 
B. 
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Theorem 2 (Algebraic structure of representations) Let V,W G V, K, S G and 
f,g & C°°(R). If two equalities 

{EV)U-UV = U'[K] + fUx, (19) 

{EW)U -UW = U'[S] + gUx, (20) 

hold, then we have a third equality 

{ElV, WJ)U - UIV, Wj = U'[T] + If, gjUx, where T = [K, S], (21) 

where {V, W] , [K, S] and |/, g\ are defined by ^TB^), ( fl^ j and (|7^j, respectively. 

According to this theorem, we can easily find that if a system Ut = K{n, u) is isospec- 
tral, i.e., At = / = 0, then the product system Ut = [K, S] for any S E B'^ can be viewed 
to be still isospectral because we have |/, g\ = [0, g] = 0, where g is the evolution law 
corresponding Ut = S{n, u) . Actually the above theorem gives a discrete zero curvature 
representation for a product system ut = [K,S], which possesses the same order matrix 



operators as ones for the original systems Ut = K(n, u) and Ut = S{n, u) (see ||2^ for the 
continuous case). Combining two theorems above can show the following. 

Corollary 1 The space defined by 

{{K, V, f) e {B\ V^ C°°(R))| U'[K] + fUx = {EV)U - UV} 
is a Lie sub-algebra of (i3^, V", C°°(R)) under the Lie product ^Tlj). 

This corollary tells us a Lie algebraic structure about zero curvature equations and will 
help us to establish Lie algebraic structures of symmetries including master symmetries. 

However, for zero curvature representations, some interesting problems remain to be 
solved. For example, assuming that two initial systems Ut = K{n,u) and Ut = S{n,u) 
have zero curvature representations possessing different order matrix operators, we want 
to know whether there exist any zero curvature representations for the product system 
Ut = [K, S] and what structures the resulting zero curvature representations possess if 
the answer is yes. It is likely to be helpful in solving this problem to use the Kronecker 



product as [g2 



III Lax operators and master symmetries 

Assume that we already have a hierarchy of isospectral integrable systems of discrete 
evolution equations of the form 

Ut = Kk = ^''Ko, $ G Kq G S^ k> 0. (22) 
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or of the form 



Ut = Kk = JGk = MGk-i, J, Me V, Gk-i eB\ k> 0. (23) 

associated with a discrete spectral problem 

Ect> = Uct>, = (0i,---,0,r. (24) 

The second form (^) more often occurs than the first form (0), although it is simpler 
to deal with the first form (p2D. Generally speaking, the operator $ above is a hereditary 
symmetry operator (see p3| for definition) determined by the spectral problem (p^ ) and 



J, M constitute a bi-Hamiltonian pair ||2^ If we choose $ = MJ~^ when J is invert- 
ible, then the form ( p3D may be changed into the form ([2^). Usually $ involves nonlocal 
operators, for example, A~^, but J, M often involves only local operators. Our examples 
are all local Hamiltonian systems. 



III.l Structures of Lax operators 

For a given X e B'^ or G G B'^ , let us introduce an operator equation of f2 G V^: 

{En{x))u - un{x) = u'l^x] - xu'[x], (25) 

in the case of (P^), or an operator equation of flj G V^: 

{EQj{G))U - UQj{G) = U'[MG] - XU'[JG], (26) 

in the case of (^). We call them the characteristic operator equations of U. The intro- 
duction of the operator equation (pSj) (or (pB])) is an important step in our manipulation. 
Obviously, we can choose flj{G) = Q{JG) when $ = MJ~^. We demand that ( pS]) (or 
(p6|)) has solutions, and VL = fl{X) (or Qj{G)) is a particular solution at X (or at G). 
Usually ( p5D (or (|26|)) has infinitely many solutions once one solution exists, because we 
can construct others fi(X) + fV for any / G C°°(R) when V E V C[A, A"^] solves 
the stationary discrete zero curvature equation {EV)U — UV = 0. The existence of solu- 
tions of {EV)U — UV = may be resulted from the existence of an isospectral hierarchy 
associated with Ecp = Ucp. 



Theorem 3 (Structure of Lax operators) Let two matrices Vq, Wo G and two vector 
fields Ko,po eBi (or po = J70, 70 G B'^) satisfy 

{EVo)U-UVo = U'[Ko], (27) 

{EWo)U - UWo = U'Ipo] + XUx. (28) 
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// we define pi, I > I, Vk, k > 1, and Wi, I > 1, as follows 

Pi = $Vo, l>l{or pi = J-ii = M-ii.u 11 eB\ 1>1), (29) 

k 

V|, = X% + J2^'"'^i^^^-l) {or f2j{G,-,)), k>l, (30) 



i=l 
I 



Wi = X'Wo + J2^'''^iP,-i) (or 1>1, (31) 

then Vk,Wi, k,l > 0, satisfy 

{EVk)U -UVk = U'[Kkl {EWi)U -UWi = U'[pi] + \'+^U^, kj>0. (32) 

Therefore for any k,l > 0, the systems of discrete evolution equations Ut = Kj. and Ut = pi 
possess the isospectral (Xt = 0) and noniso spectral (Xt = X''^^) discrete zero curvature 
representations 

Ut = {EVk)U - UVk, Ut = iEWi)U - UWu 

respectively. 

The theorem shows that the Lax operators associated with two hierarchies of interest- 
ing vector fields can be constructed simply by a unified form. Its proof is left to Appendix 
C. We are successful, thanks to introducing a characteristic operator equation. The dif- 
ficulty is now transferred to seeking a solution to the characteristic operator equation. 
However this can automatically be solved on basis of the structure of Lax operators of 
isospectral hierarchies, which will be seen in the next subsection. 

III. 2 A method for constructing master symmetries 

Now we focus our attention on the construction problem of master symmetries. Theo- 
rem R| already shows the structure of Lax operators associated with the isospectral and 



nonisospectral hierarchies (refer to for the continuous case). When an isospectral hier- 
archy ( p2D (or (p3D) is known, the theorem also provides us with a method to construct a 
nonisospectral hierarchy associated with the discrete spectral problem ( ^4]) by solving an 
initial discrete zero curvature equation (^) and solving a characteristic operator equation 
(B5I) (or m). 



However, a solution to (|25| ) (or (0)) may easily be generated by observing the resulting 
Lax operators. In fact, we have 

^{Kk) (or nj{Gu)) = Vk+i - AV^fc. (33) 
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This may be checked, say, for the case of p^), as follows 

by using (|30|). Now by the first equality of (0), we may compute the following 

{En{Kk))u - un{Kk) 

= {EVk+i - \EVk)U - U{Vk+i - XVk) 

= {{EVk+i)U -UVk+i)-\{{EVu)U -UVu) 

= U'[Kk^^] - \U'[Kf,] = U'[^Ku] - Af/'iiTfc], 

for example, for the case of (|^). Therefore we see that a possible solution VL{X) to (|25|) 
(or VLj{G) to ([26|) ) may be generated by replacing the element Kk (or Gk) in the equahty 
(H) with X (or G). 

The Lax operator matrices Vk+i and Vk are known, when the isospectral hierarchy 
has already been found. Thus we don't have to directly solve the characteristic operator 
equations and then the whole process of construction of the nonisospectral hierarchy 
becomes an easy task, finding po, to satisfy (^) and computing Vk+i — AVfc to find a 
solution to (^ (or (^)). 

The nonisospectral hierarchy (^) is exactly the master symmetries that we need to 
find. The reasons are that the product systems between the isospectral hierarchy and the 
nonisospectral hierarchy are still isospectral by Theorem ^ or as we said before in Section 
II, and that usually all systems of the isospectral hierarchy commute with each other. 
Therefore it is because there exists a nonisospectral hierarchy that there exist master 
symmetries for isospectral systems of discrete evolution equations derived from a given 
discrete spectral problem. 

In the next section, we shall in detail illustrate our construction process by three 
concrete examples and establish the corresponding centerless Virasoro symmetry algebras. 



IV Applications 

We illustrate only by three examples how to apply the method in the last section to 
construct master symmetries for various lattice hierarchies. 

To make the process clearer, we introduce a conception for a given discrete spectral 
problem Ecf) = U(f), which has an injective Gateaux derivative U'. That is a uniqueness 



property similar to the one in the continuous case [p7[] : 



if {EV)U -UV = U'[K], V eV^^ C[A, A'^], K G B% and V\^=o = 0, then V = 
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and further = by the injective property of U'. It means that if an isospectral (A^ = 0) 
Lax operator V equals zero at m = 0, then so does V itself. Actually, this property cor- 
responds to the uniqueness of an integrable hierarchy associated with a spectral problem 
E(j) = U(j). That is to say, when initial conditions and constants of inverse difference op- 
erators are fixed (for example, as in (|^) and (^), the associated isospectral hierarchy is 
uniquely determined. Most of discrete spectral problems share the uniqueness property. 
The following three spectral problems are exactly examples which share such a property. 

IV. 1 The Volterra lattice hierarchy 



Let us first consider the following discrete spectral problem [15 



The corresponding isospectral integrable lattice hierarchy reads as 

Ut = K, = $^^0 = w(4'+i - 4+1 ), Ko = «(«(-!) - ^(1)), A; > 0. (35) 
/ ^ (1) \ 

Here the matrix V = J2i>o * *+i solves the stationary discrete zero curvature 

- \ Ci -tti J 

equation {EV)U — UV = 0, where we choose the initial conditions 

1 

a-o = 2' ^ ^'^ ^ ^ 
and the hereditary operator $ is given by 

$ = m(1 + E-'^){-u^^^E'^ + u){E- (36) 

where {E — 1)~^ is determined by (|^). It is worth pointing out that each system in ( p5D is 
local and polynomially dependent on u, although the hereditary operator $ has nonlocal 
and non-polynomially dependent features. 



The first discrete evolution equation is the Volterra lattice equation 

{u{n))t = u{n){u{n — 1) — u{n + 1)), 

which is significantly generalized by Bogoyavlensky [^]. The associated Lax operators are 
as follows 

V, = (A'=+V)>i + I ], k>0, (37) 
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where {P)>i denotes the selection of the terms with degrees of A no less than 1. In 
particular, the initial isospectral Lax operator reads as 

}:X — u Xu 

1 -iA-.<-.. 1- '^^ 



The result till here can be obtained from (|34D by using a powerful method in [^ . 

We easily obtain the corresponding quantities in the nonisospectral (A* = A) initial 
discrete zero curvature equation (]28|): 

p, = u,Wo=(^^^ j , (39) 

and a solution to the characteristic operator equation (^) by (^): 

m) = \ ' J " J 1, (40) 



where Qij{X), i,j = 1,2, are given by 



fin(X) = {E- 1)-1(-m«E2 + u)iE- l)-'u-^X 
n^iX) = \uE{E - i)-^u-^x 

VL2i{X) = {E-l)-^u-'^X 

fi22(X) = -X{E - l)-^u-^X + E-\E - 1)-1(-m«E2 + u){E- l)-\~^X. 

Now by Theorem ^ we obtain a hierarchy of nonisospectral discrete evolution equations 
Ut = Pi = $'pO) ^ > 0, associated with the spectral problem (^). 

Let us now consider how to compute the corresponding symmetry algebra. The idea 
below can be applied to other cases. We first make the following computation at m = 0: 

Kk\u=Q = 0, pi\u=o = $VoU=o = 0, k,l>0, 
Vk\u=o = X'(^ _iA ) ' ^ - °' 



2 



Vkx\u=o = kX'-' I ) ^ ( -i ' ' ^ - °' 

Wix\u-o = ^A'-^ I ^ ° I + (1 - 610) I ° , ° 1 , / > 0, 

' -2 J \{l-l)X'-^[n] -lX'-'[n] - ' 
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where Vk, Wi, k,l > 0, are given as in Theorem |^ and 5/o represents the Kronecker symbol. 
While computing py/|u=o, we need to note that fi(po)|n=o 7^ 0, but fl{pi)\u=o = 0, Z > 1. 
The other two examples below have a similar character, too. Now we can find by the 
definition (0) of the product of two Lax operators that 

lVk,Vi]\u=o = 0, k,l>0, 

lVk,Wi]\^=o=ik+l)Vk+i\u=o, kj>0, (41) 
lWk,Wi]\u=o = {k-l)Wk+i\u=o, kj>0. 
For example, we can compute that 

lVk,Wi]\^=o = [Vk\u=o,Wi\u=o] + X'^'Vkx\u=o 

\ / i 




(A; + 1)A 



Because [Vk, Vi], [Vk, Wi] - (A; + 1)\4+«, {Wk, Wij -{k- l)Wk+i, k,l>0, are all isospectral 
{\t = 0) Lax operators belonging to ® C[X, X~^] by Theorem ||, based upon (|l|) we 
obtain a Lax operator algebra by the uniqueness property of the spectral problem (|3^) 

lVk,Vij = 0, k,l>0, 

lVk,Wij = {k+l)Vk+i, k,l>0, (42) 
lWk,Wil = {k~l)Wk+i, kj>0. 

Further, due to the injective property of U', we finally obtain a vector field algebra of the 
isospectral hierarchy and the nonisospectral hierarchy 

[Kk,Ki] = 0, kj>0, 

[Kk,pi] = ik + l)Kk+i, k,l>0, (43) 

[Pk, Pi] = {k- l)pk+u k,l>0. 

This implies that pi, / > 0, are all master symmetries of each lattice equation Ut = Kj^^ 
in the isospectral hierarchy, and the symmetries 

Kk, k>0, and r/^°^ = t[Kk„ pi] + pi, I > 0, 

constitute a symmetry algebra of Virasoro type possessing the same commutator relations 

as m. 
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IV. 2 The Toda lattice hierarchy 



Let us second consider the discrete spectral problem |5D 

E(j) = U(f), U = 



1 

—V X — p 



u 



P 



(44) 



The corresponding isospectral integrable Toda lattice hierarchy |31] reads as 



V — V 



(1) 



k>0. (45) 



Here V = 



i>0 



0-i bi _. 

(^-) I A^* solves {EV)U — UV = 0, in which we choose 



—vb] —Qi 



1 



oo = ^0 = 0, ai = 0, bi = -1, 



and the hereditary operator $ is determined by 
$ = 



p {v^^'^E^ ~v){E -l)-^v-^ 

v{E-^ + l) v{pE - p^-^^){E - l)-^v-^ 



(46) 



The first system of discrete evolution equations is the Toda lattice 

{p{n))t = v{n) - v{n + 1), 
{v{n))t = v{n){p{n) - p{n - 1)), 

up to a transform of dependent variables. The lattice hierarchy above has a local tri- 
Hamiltonian structure 

SHk+2 SHk+i .r^Hk 

ut = Kk = J— = M— = N^—, fc > 0, 

ou ou ou 

where the Hamiltonian operators J, M, N and the conserved quantities if^ defined by 



J 





v{E-^ - 1] 



{1-E)v 




AT = M^^ = -<I>M 

_/ p{vE~^ - Ev) + {vE-^ - Ev)p p'^{l-E)v + {vE-^ -Ev){l + E)v 
~ \ v{E-^ + l){vE-^ - Ev) + v{E-^ - l)p'^ 2v{E-^p - pE)v 

Ho=p + ^Inw, Hk = /c > 1, 

2 k 
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where denotes the conjugate operator of $. Note that this tri-Hamiltonian structure 
may be estabhshed through a trace identity ||3^. The corresponding Lax operators read 



as 

V. = iX'^'V), + ( iyk> 0, (47) 

where the subscript + denotes to select the non-negative part. Hence in particular 



It is easy to find the corresponding quantities in the nonisospectral {Xt = X) initial 
discrete zero curvature equation (]28|): 

where [n] is the multiplication operator defined by and a solution to the characteristic 
operator equation by (P^): 

where flij{X), i,j = 1,2, are given by 

nn(X) = E-\E - l)-'Xi + - X){E - l)-^-^X2, 

nu{X) = {E-l)-'v-'X2, 

n2i{X) = vE{E - i)-\-^X2, 
n22{x) = {E-i)-'x,. 

In this way, we obtain a hierarchy of nonisospectral systems of discrete evolution equations 
Pi = ^''Po, I > 0, associated with the spectral problem (|^). 

In order to construct a vector field algebra, we make a similar computation at u = 0: 



Kk M=o 


— 0, pl\u=0 






^fc M=0 = 























A; > 0, 



-2 AM 2M 
\n\ I 



U -1 \ , / i 



-iA; \^ -\ 
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Now we can find through the product definition of [■, ■] in ( ]T8| ) that 

lV^fc,V^lk=o = 0, k,l>0, 

lVk,Wi]\u=o = (k + l)Vk+i\u=o, k,l>0, (51) 
lWk,Wi}\u=o = {k-l)Wk+i\u=o, k,l>0. 

A similar argument yields a Lax operator algebra by the uniqueness property of the 
spectral problem ( Pf ) 

lVk,Vij = 0, k,l>0, 

lVk,Wi} = {k + l)Vk+i, kj>0, (52) 
lWk,Wil = {k-l)Wk+i, k,l>0. 

And then because of the injective property of U', we obtain a semi-product Lie algebra 
of the isospectral hierarchy and the nonisospectral hierarchy 

[Kk,Ki]=0, k,l>0, 

[Kk,pi] = {k + l)Kk+i, k,l>0, (53) 
[Pk, Pi] = {k- l)pk+i, k,l>0, 

which gives rise to a symmetry algebra of Virasoro type for the isospectral Toda hierarchy 



IV. 3 A sub-KP lattice hierarchy 



Let us finally consider the discrete spectral problem |^| 



E^ = U(p, U 



/ 1 0\ 

b — X a 1 
V c 0/ 



u 



b 

V c ; 



V 03 y 



(54) 



which is equivalent to {—E'^ + b + aE + E ^c)0i = A0i, a sub-KP discrete spectral problem 
^M. The corresponding isospectral integrable lattice hierarchy reads as 



ut = Kk = JGk = MGk-i, k>0, 
where a Hamiltonian pair J, M and Go, Gi are defined by 
( E-E-^ \ 



(55) 



J 



\ 



(E-i - l)c 

-c{E-l) 
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/ Eb- bE-^ + aA+A-iA_a EcE - E'^c 
M = cE-E-^cE-^ E-^ac-acE -bA^c 

-cA+b c[A+A-iA_ - A_ - A+]c / 
c{Eb + b) 

G-i = I 1 |,Go=| b |,Gi= b^ + ac + E-^ac 



cA+ - 


A-iA_a 




fc\ ( 


j , Go = 


b , Gi = 




\a 1 \ 



where A+, A_ are the difference operators: A_| 
system of discrete evolution equations is 



E-h A_ 



1 — E ^ . The first nonlinear 



{a{n))t = c{n + 1) - c{n - 1), 
< {b{n))t = a{n — l)c{n — 1) — a{n)c{n), 
{c{n))t = c{n){b{n) — b{n + 1)). 

We easily find the corresponding quantities in and 



( {E~E~')c \ 

{E-^ - l)ac 
[ c{l-E)b ) 



( 



Vo 



1 

c 

y -E~^ac E-^c X-b J 



Po = Jlo = M-f_i = J 
( 1 



.-1 A-1/ 



M 



\n\ 







\ -c-^A-_^b j 









\ 



KM + i) 







V KM + 2) ; 

We can also obtain a solution to the characteristic operator equation (123) by (P3|): 













( ^^'^ ] 




1^21 (G) 




^^23(G') 


, G = 


G{2) 








^zz{G) ) 




\ G(3) ) 



where flij{G), i,j = 1,2,3, are determined by 

Qu{G) = -{E' + E)-\cG(^s^ + EaG^.j) 
^uiG) = E 

^13 (C) = ^(2) 

fi2i(G') = cEG(2) + (6-A)G(i) 



^^22(G') = -{E + l)-\cG^3) + ^aG(i)) + aG 

^23 (C) = 



(1) 



(56) 



(57) 
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Q^siG) = -E{E + l)-\cG(3) + EaG(^^) + A+aG^) - {b - A)G'(2). 

By Theorem ^ we get a hierarchy of nonisospectral systems of discrete evolution equations 
Ut = Pi = $'pO) ^ > 0, associated with the spectral problem (^). 

In order to generate a vector field algebra of the isospectral hierarchy and the non- 
isospectral hierarchy, we need the following quantities, which may be directly worked 
out: 

Kk\u=o = 0, poU=o = -''7oU=o = 0, pi\u=o = Jli\u=o = M-fi_i\u=o = 0, /c > 0, / > 1, 
/ 1 \ 



VfcL=o — A 





V A y 







f i 
















/ 




-IM ^ 








= A' 





i(M + i) 







+ (1 


- Sio)X'- 


-1 


















lo 





KM - 


f 2) ) 








^ 




-f AH ) 












/ 1 ^ 




/ \ 


















VfcA M=0 


= kX 


fc-1 





+ A*^ 



























y A yl 




[o l) 
































\ 











-w- 


- 1)A'-2H 


\ 


Wix\u=0 


= IX 


l-l 


i([n 


] + l) 





+ (1- 























V 





KM 


+ 2) 


/ 






lo 





-\ix^-M 


/ 



Now we easily find according to the product definition of [■, ■] that 

[Vfe,^lk=o = 0, fc,/>0, 
{Vk,Wi\\u=^ = {k + l)Vk+i\u=^, k,l>0, 
lWk,WijU=o = {k-l)Wk+i\u=o, k,l>0. 
The same deduction leads to a Lax operator algebra 

lVk,Vi} = 0, k,l>0, 
lVk,Wij = {k + l)Vk+i, fc,/>0, 

m,wi} = {k-i)Wk+i, k,i>o. 

and further a vector field algebra 

[Kk,Ki] = 0, k,l>0, 
[Kk,pi] = ik + l)Kk+i, k,l>0, , 
[Pk, Pi] = {k- l)pk+i, k,l>0. 



(5^ 



(59) 
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which may generate a master symmetry algebra possessing the same algebraic structure 
as (|9l). 



V Conclusion and remarks 



We have established an algebraic structure related to discrete zero curvature equations 
and further introduced a simple but systematic approach for constructing master sym- 
metries of first degree for isospectral lattice hierarchies associated with discrete spectral 
problems. The resulting algebraic structures also leads to an explanation of why there ex- 
ist master symmetries of first degree. Some complicated calculation in our construction is 
saved by using a characteristic operator equation (|25| ) (or (|^)) and a uniqueness property 
of discrete spectral problems. The crucial step is the construction of the corresponding 
nonisospectral lattice hierarchies, which can be found by solving an initial nonisospectral 
discrete zero curvature equation. Three lattice hierarchies are showed as illustrative ex- 
amples and the corresponding master symmetry algebras of centerless Virasoro type are 



exhibited. Some of the results in this paper have been reported at SIDE 11, UK [35 



It is worth noting that three examples described in the last section possess the same 
commutator relations between their isospectral and nonisospectral vector fields. In gen- 
eral, we have [Kk, pi] = {k + 'y)Kf:+i, 7 = const., but the other two equalities of the whole 
Virasoro algebra don't change. This is also a common phenomenon for continuous inte- 
grable hierarchies \^\ ^7^. Furthermore we may add a nonisospectral master symmetry 
with At = 1 to the whole Virasoro symmetry algebra but this often requires additional 
checking. For example, a nonisospectral master symmetry with A* = 1 of the sub-KP lat- 
tice hierarchy (|55|) is p_i = J7-1 = (0, 1, 0)"^. On the other hand, similar to the theory in 
|37| , we may also choose an operator solution Q{X) (or Qj{G)) satisfying Q{X)\x=o = 
(or Qj{G)\g=o = 0) (all three examples in the last section have this property) and then 
we only need to compute [Vq, VFo]|u=o so as to give Lax operator algebras at m = and 
finally give Lax operator algebras generally. 

In our discussion, in fact, we haven't used the hereditary property of the recursion 
operator $ (or the bi-Hamiltonian property of J and M), while we construct Virasoro 
symmetry algebras for integrable lattice hierarchies, and thus it can also be applied to 
lattice hierarchies which possess non-hereditary recursion operators. The advantage of 
our scheme is to fully utilize discrete zero curvature equations so that the whole process 
to generate master symmetries of first degree becomes an easy task. There were also 
an algorithm implemented in MuPAD and other direct tricks [|I^ [|T3| to 
compute master symmetries of first degree for systems of discrete evolution equations. 
However our theory focuses on seeking an answer to the existence and structure problem 
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of master symmetries of first degree. 

We should mention that there exists a large variety of other theories or methods to dis- 
cuss integrable properties of systems of nonlinear discrete equations, which include Hamil- 
tonian theory |^0|] Backlund-Darboux transformation |^ R-matrix method 

symmetry reduction |^ etc. Moreover we can consider the time discretization prob- 



lem and periodic initial and boundary value problems of time discretizations 
for symmetry flows of systems of discrete evolution equations. The resulting difference 
equations and mappings should be useful in discussing the integrability of the underly- 
ing systems of discrete evolution equations themselves. We are also curious about the 
following natural problem: Are there any higher degree master symmetries for systems 
of discrete evolution equations which don't depend explicitly on the evolution variable? 
If the answer is yes, can we establish any relations between those higher degree master 
symmetries and discrete zero curvature equations as we did for the first degree master 
symmetries? 



Appendix A: Proof of Theorem [1] 

Let {Ki,Vi,fi) e V^C~(R)), 1 < z < 3. Because the bilinearity and the skew- 
symmetry of the product (|T^ are self-evident and we already know that the products 
defined by ( ]T^ and (|T3]) are Lie products, we only need to prove the following Jacobi 
identity: 

[[^i,V^2l,V^3l+cycle(l,2,3) = 0. (A.l) 
Let us first compute by ([T8| ) that 

= {v;[K,]y[Ks] - {v^[K,]y[Ks] + [v,, v^lli^s] + f2{v^xy[Ks] - fi{v,,y[Ks] 

-Vi[[K,, K2]] + [V;[K2], V,] - [V^[K,], Vs] + [[^1, V2],Vs] + ^[I^ia, Vs] - h[V2x, V^] 
+MV;[K2])x - f3mK,])x + fs[V,, V2]x + f2xf3Vix + /2/3V1AA 

-/ia/3V^2A - /l/3^2AA " I/l, /2F3A. (A.2) 

We need to use the following fundamental equalities 

iVxy[K] = iV'[K])x, V G V, K e B\ 

[V, W]x = [Vx, W] + [V, Wxl V,We V^ 

[V^, W][K] = [V'[K],W] + [V, W'[K]], V,WeV, K e B\ 

v'[T] = iy'[K]y[s] - {y'[s]y[Ki t = [k,si v e v^ K,se b\ 
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which may be shown by a direct computation and the last equality of which is a similar 
result as in [^. Now we can go on to compute that 

= {vi[K,]nKs] - mK.mKs] - mK,]y[K2] + mK2]y[K,], 

A,''' := [V^,V2nKs] + [Vl[K,],V,]-[V^[K^],Vs] 

= [v;[K,],v2] - [v^iK^i Vi] + [v;[K,],v3] - Mm,V3], 

Ac''' := /2(V^lA)'[i^3]-/l(V^2A)'[i^3]+/3(K[^2])A-/3(V^2[^l])A 

= f2{V,xy[Ks] - fliV^xYlKs] + /3(V^lA)'[i^2]) - /3(V^2A)'[i^l], 

A/'' := /2[1^1A,V^3]-/i[V^2A,V^3] + /3[V^1,V^2]a 

= /2[1^1A, Vs] - h[V2X, V3] + /3[FlA, V2] - f3[V2X, Vi], 

Ae''' := /2a/3V^1A + /2/3V^lAA - /ia/3V^2A " /l/3l^2AA " I/l, /2F3A, 

= /2A/3V1A + 72/3 VlAA — /1A/3V2A — /1/3V2AA — /1A/2V3A + /i/2aV3A. 

A direct check can result in that 

A*"'^^'^ + cycle(l, 2, 3) = 0, where * = a,b, c, d or e. 
Noting ( |A.2| ), it follows therefore that 



II^i,^2l,^3l + cycle(l,2,3) 
= A^''' + A,''' + Ac''' + A/'' + Ac''' + [[V„ V2], Vs] + cycle(l, 2, 3) = 0, 

which is exactly the Jacobi identity (|A. 1|) and thus completes the proof. 



Appendix B: Proof of Theorem |2] 

The proof is an application of the equalities (|TUp and ( PDj ) and the third equality 

{U'[K]nS] - {U'[S]y[K] = U'[T], T = [K, S], (B.l) 

which has been mentioned in the proof of the first theorem. We observe that 

(Eqn. ® - (Eqn. mYiK] + 9{Eqn. (p])), - /(Eqn. 

The resulting equality reads as 

mK]y[S]-iU'[S]y[K] + lf,g]U, 
= {EV'[S])U + {EV)U'[S] - U'[S]V - UV'[S] 

-{EW'[K])U - {EW)U'[K] + U'[K]W + UW'[K] 
+g{EVx)U + g{EV)Ux - gU^V - gUV^ 

-f{EWx)U - f{EW)Ux + fUxW + fUWx. (B.2) 
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On the other hand, we have immediately 

{ElV,W])U-UlV,Wj 
= {EV'[S])U - {EW'[K])U + {EV){EW)U - {EW){EV)U 
+g{EVx)U - f{EWx)U - UV'[S] + UW'[K] 

-UVW + UWV - gUVx + fUWx- (B.3) 
It follows therefore from (|RlD , O and that 

{ElV, Wl)U - UIV, Wj - U'[T] - If, gjUx 
= (ElV, Wj)U - UIV, Wj - mmnS] + mS]y[K] - U, gjUx 
= {EV){{EW)U - V'[S] - gUx} - {EW){{EV)U - U'[K] - fUx] 

-UVW + UWV + gUxV - fUxW + U'[S]V - U'[K]W 
= {EV)UW - {EW)UV - UVW + UWV + gUxV - fUxW + U'[S]V - U'[K]W 
= {{EV)U -UV - fUx - U'[K]}W - {{EW)U - UW - gUx - U'[S]}V 
= 0, 

which is what we need to prove. 

Appendix C: Proof of Theorem g 

We prove two equalities in (p2D. The rest is obvious. We compute that 
iEVk)U-UVk 

k 

= X'[iEVo)U - UVo] + J2^'~'[E^iK^-l)U - UQ{K,^i)] 

i=l 

= X'^U'lKo] + ^ X''-'{U'[<^>K,.^] - XU'[K,_^]} 

i=l 

= X'U'lKo] + E >^'-'{U'm - XU'[K,^^]} 

i=l 

= U'[Kk], k>l; 
{EWi)U - UWi 

I 

= X'[{EWo)U - UWo] + E x'-^[En{pj^,)u - un{pJ^^)] 

i=i 

= A'{f/'[po] + XUx} + E V-^{f/'[<l>p,_i] - Af/'[p,_i]} 
= X'{U'[po] + m} + E A'-^{f/'[p,] - XU'[p,.,]} 
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Note that we have used the characteristic operator equation (P3| ) but the situation in the 
case of (^) is conpletely similar. The proof is therefore finished. 
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